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ABSTRACT 


We have derived the radial distribution of gamma-ray emissivlty In the 
Galaxy from flux longitude profiles, using both the final SAS-2 results and 
the recently corrected COS-B results and analyzing the northern and southern 
galactic regions separately. We have then made use of the recent CO surveys of 
the southern hemisphere, In conjunction with the northern hemisphere CO data, 
to derive the radial distribution of cosmic rays on both sides of the galactic 
plane. We have found that. In addition to the "5 kpc ring" of enhanced 
emission, there Is evidence from the asymmetry In the radial distributions for 
spiral features which are consistent with those derived from the distribution 
of bright HI I regions. We find positive evidence for a strong Increase In the 
cosmic ray flux In the Inner Galaxy, particularly In the 5 kpc region. In both 
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I. INTRODUCTION 

It has long been recognized that Y-ray astronorty provides a powerful tool 
for studying the origin and distribution of the primary cosmic radiation 
(Ginzburg and Syrovatsky 1964; Fazio 1967) and also for revealing the 
structure of the Galaxy (Stecker 1971, 1977, 1978). An important feature of 
large scale galactic structure was obtained quite early from analysis of the 
preliminary SAS-2 data when it was deduced that the Y-ray emissivity of the 
Gala>;y has a large maximum in a ring at about 5 kpc from the galactic center 
(Stecker, et al. 1974) and this feature was then quickly related to a newly 
discovered 5 kpc ring of giant molecular clouds in the Galaxy (Solomon and 
Stecker 1974, Scoville and Solomon 1975, Stecker et al. 1975). 

Because the Galaxy is transparent to Y-rays, Y-ray surveys can be used in 
conjunction with galactic surveys at other wavelengths to provide a "synoptic" 
approach to the general problem of galactic structure (Stecker 1981). There 
have been basically two distinct approaches for theoretical research on this 
problem. One approach, which can be called the "modeling approach", is to 
adopt a model for galactic spiral structure based on observations at more 
traditional astronomical wavelengths, to then make an assumption relating the 
cosmic ray and gas distributions, and to check for consistency between the y- 
ray distribution implied by the model and the actual observations. This 
approach has been used extensively by Fichtel, Kniffen and coworkers (e.g., 
Fichtel and Kniffen 1984). 

The other approach is to start with the Y-ray observations themselves, 
and to use them together with information gathered at other wavelengths to 
deduce new information about galactic structure and the galactic cosmic ray 
distribution by employing geometrical unfolding and analysis techniques (Puget 
and Stecker 1974, Solomon and Stecker 1974, Stecker 1975, Strong 1975, Stecker 
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Caraveo and Paul 1979). In this approach, a key element has been the 
understanding that molecular hydrogen plays an important role in the problem 
(Stecker 1969, Stecher and Stecker 1970), Because of this fact, applications 
of the unfolding approach have been limited by the lack of information on the 
distribution of molecular hydrogen on the southern hemisphere side of the 
galactic disk. This lack has now been remedied as the results of recent 
southern hemisphere survey are now becoming available (Sanders, Solomon and 
Scoville 1984; Robinson et. al. 1984). In this paper, we make use of these 
new results, as well as the final SAS~2 Y-ray survey data, and the COS-B Y-ray 
survey results, newly corrected for intrinsic background, in order to 
reexamine the galactic structure and galactic cosmic-ray gradient problems. 

II. GAMMA-RAY FLUX UNFOLDING METHOD 

Observations of the galactic gamma-ray emission in the form of flux- 
longitude profiles give the total integrated emission along a line of sight. 
Several methods have been used to unfold the flux-longitude profiles and 
derive the gamma-ray emissivity as a function of galactic coordinates, Q(R,0), 
making a necessary assumption about the synnnetry of Q{R,0). Strong (1975), 
Strong and Worrall (1976) and Caraveo and Paul (1979) used a numerical matrix 
inversion to unfold the SAS-2 data, assuming either cylindrically symmetric or 
spiral arm models for Q. Puget and Stecker (1974) derived an analytic 
expression for the unfolded radial gamma-ray emissivity distribution, assuming 
cylindrical symmetry. They then proceeded to a numerical analysis using the 
earliest SAS-2 data. We will use this latter method to unfold the most 
recently published results from both the SAS-2 and COS-B satellite detector 
experiments. 

We assume cylindrical symmetry in each half of the galactic plane so that 


Q(R,0) = Q(R) Is a function only of galactocentric radius R, and is 
independent of distance 2 above the galactic plane up to a height, h. The 
gamma-ray flux, in photons cm"^ rad"^ can be written. 


I (A) 


«0 


K .. I (h/l!0)cotl> 
>0 ^ 


dp Q(r) 


( 1 ) 


where Rq is the galactocentric distance of the Sun, p is heliocentric 
distance, r « R/Rq and b is galactic latitude. We divide I(t) into flux 
contributions from the inner (r < r^) and outer (r > r^) parts of the Galaxy, 
I =» + Iq, and assume that the emissivity in the outer Galajy, Q^(R)» is 

constant for rQ < r < r^,. Using Eq (1) and p = cosJi + (r^ - sin^A)!/^^ 
find the following expression for the outer gala;y flux component 



where, (2) 

P|ji = cosA - (r^ - sin^Jl)^'^^ 

p“ = cos^i ± (r^ - sin^A)^'^^ 


Using this expression, we can determine the part of the observed flux 
from emission at galactic radii less than rg and invert this flux to find the 
inner Galaxy emissivity Q^(r). Using Eq (1), the inner galaxy flux may be 
written. 
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dr 
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(3) 


Puget and Stecker (1974) have shown that this equation can be Inverted with 
the use of Laplace transforms to give the following solution for Q^(r), 

2 




(4) 


Where n = sln^A. The cos^ In the denominator of the Integrand accentuates 
local fluctuations In tho flux (i- 90“ and a 270") which can overpower the 

flux derivatives In more distant regions. Since we are mainly Interested here 
In the large-scale distribution of gamma-ray emissivlty In the Galaxy, we can 
eliminate these fluctuations In the unfolding due to local sources of emission 
by taking rg = 0.85, which corresponds to longitudes 300" < a < 60*. 

We have used this technique to unfold both SAS-2 and COS-B longitude 
data. The SAS-2 data for energies > 100 MeV from Hartman et. al (1979) Is 
Integrated over latitudes |b| <_ 10" and has a longitude resolution of 2.5". 
Mayer-Hasselwander (1983) has presented a COS-B longitude profile for energies 
> 100 MeV with a longitude resolution of 2.5" for direct comparison with the 
SAS-2 profile. With a COS-B background subtraction of 8 X 10“ ^ ph cm"^ s"^ 
sr”^ above 70 MeV, the two data sets were found to be In agreement within the 
statistical errors at most longitudes. It was also found that the COS-B flux 
must be multiplied by a factor of 1.17 to bring the general flux levels of the 
two data sets Into agreement. In addition to this COS-B data set, we have 
unfolded the COS-B data In the energy range 0.3 - 5 GeV from Mayer- 
Hasselwander et. al. (1982), which were also averaged over (bl _< 10" but have 
a longitude resolution of 1". 

An outer galaxy flux contribution Is determined from Eq (2), assuming a 
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constant value » 1.1 X 10"25 cm“3 s“i (E > 100 MeV) and Qq = 0.51 X 10‘2S 
cm"^ s"M.3 < E < 5 GeV) and constant scale height h = 150 pc between Tq = 0 
.85 and r,^ « 1.4. The value of Qq = qQ <n^> was derived from a local gamma- 
ray production rate qQ = 1.9 X 10~2S g-i fo^ g > iqo MeV (Stecker 1977, 

Fichtel and Kniffen 1984) and qg = 8.8 X 10"26 g-l fop 0.3 < E < 5 GeV, taking 
a local average H-atom density of <n[^> = 0.6 cm“3. This outer galaxy 
emissivlty gives approximately the flux levels observed In the anticenter 
directions. Our adopted values for qo agree with those recently derived by 
Sloemen et al. (1984) to within 15%. They obtain qg (>100 MeV) = 2.2 x lO-^s 
s"^ (extrapolating from a lower energy of 70 MeV) and q^ (0.3-5 GeV) « 7.4 x 
10“25 s"^ 

Before unfolding, the flux data points were first averaged over a 
longitude range of 7.5“ in the case of the > 100 MeV data and 7“ In the case 
of the high energy COS-B data In order to smooth out fluctuations. The outer 
gala;^ flux contribution, determined by the method described above, was 
subtracted from the total flux at each longitude to obtain an Inner gala^^y 
flux contribution. The derivatives d/dn(-I^/cosA) In Eq (4) were then 
evaluated from a cubic spline fit to the averaged flux data points In the 
longitude range 300" < A < 60" from which the outer gala)^ flux as determined 
above was subtracted. 

Errors In the calculated emissivltles were estimated by assuming that the 
standard deviations of the derivatives, which would be difficult to determine 
exactly from the spline fits which correlate all the points and their errors, 
are just equal to those of the fluxes at those points. Propagating these 
errors In Eq (4) then gives error values for each of the derived 
emissivltles. In order to confirm that reasonable estimates have been made 
for the errors, we have obtained the expected statistical spread In the 
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ein-|jslvlty profiles for several cases by computing different realizations of 
the profiles from a Monte Carlo sampling of the observed flux errors. The 
variation In the profiles was found to agree quite well with the error 
propagation method, which gives the actual error values and thus enables us to 
plot error bars on the emissivlty profiles. We note, however, that because of 
the correlations Introduced by the spline fits, the errors In adjacent Q(r) 
values are strongly correlated and thus the error bars do not represent 
Independent errors. 


III. RADIAL EMISSIVITY DISTRIBUTIONS 
Results from unfolding the three flux-longitude profiles de'scribed In 
Section II are shown In Figures 1-3. The radial distributions of gamma-ray 
emissivlty (times scale height) for the SAS-2 and COS-B data at energies > 100 
MeV are plotted In Fig. 1. The appearance of negative emissivltles In several 
regions, the only significant one being the 1.5 - 3.5 kpc region In the South, 
results from a breakdown In the assumption of cylindrical symmetry. 

Strong or local points sources and spiral structure could be responsible for 
such departures from cylindrical symmetry. Both of these features are 
probably present In the data. There Is general agreement In the shapes of the 
COS-B and SAS-2 emissivlty distributions, the dominant features being a peak 
between 5 and 6 kpc In the North and a peak between 4 and 5 kpc In the South, 
which seem to describe an asymmetric ring of emission. The emission appears 
to be most Intense on the Inner rim. There Is also a maximum emissivlty near 
the Galactic center and a secondary peak of emission around 7.5 kpc in the 
South, which Is more pronounced In the COS-B data. The position and magnitude 
of the peak near the galactic center are uncorrected for significant 
contributions from local emission near Ji=0" In the latitude range lbl<10“. 
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(For a more precise treatment of galactic center emission, see Blitz, Bloemen 
and Hermsen, 1984). 

A general construction of the gal'v:tic emission map is shown in Fig. 3. 
The COS-B emissivity levels in the 5 kpc region are significantly less than 
those of SAS-2 relative to the same local emissivity. There are several 
factors which influence the peak-to-local emissivity ratios in the unfolding 
process. A higher outer galajo' emissivity level would decrease the peak-to- 
local ratio, but it would do so by the same factor in both sets of data. If a 
higher background level were subtracted from the COS-B longitude flux prior to 
unfolding, the peak to local emissivity would increase, bringing it into 
better agreement with the results for the SAS-2 flux, which had essentially no 
background. Figure 2 shows the radial emissivity distribution obtained from 
the COS-B flux after first subtracting a small additional 0.2 X 10" cm"^ s"^ 
sr"^ intrinsic instrumental background correction (using 0.15 x lO"** produced 
similar results). Recent reanalysis of the COS-B detector background has 
indicated that such a higher background correction level of this order may be 
justified, although the precise amount of this correction is still uncertain 
(Mayer-Hasselwander, private communication). Within the errors, this 
emissivity distribution now agrees with the SAS-2 emissivity, but the peak-to- 
local ratios are still lower in both the North and South. 

Our results on the COS-B radial distribution of > 100 MeV gamma-rays 
(without additional background subtraction) give an emissivity in the 5 kpc 
ring which is 4-5 times the local value. This is significantly larger than 
the ratio of 2-3 obtained by Mayer-Hasselwander (1983) from an unfolding of 
the same data set. The discrepancy may be explained by the fact that our 
assumed value of <nn> “ 0*6 cm"^ (Burton 1976), is lower than his assumed 
value of 1 cm"^. Our resulting value for the local emissivity is about half 
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as much, producing a higher contrast to the emisstvity in the 5 kpc region. 

Ir. addition to changing the normalization of the radial distribution, a lower 
local emissivity gives a smaller flux contribution from the outer galaxy and 
therefore a larger remaining flux contribution from the inner galaxy. We were 
in fact able to reproduce the results of Mayer-Hasselwander's unfolding by 
taking our local emissivity to be twice the value given in Section II. 

Figure 4 shows the results of unfolding the COS-B flux data in the 
highest energy range, 0.3 < E < 5 GeV. The emissivity profile looks much the 
same as at lower energies, however the errors here are somewhat larger. The 
level of negative emissivity around 2 kpc in the South is not as large as at 
lower energies, even though the flux data has been averaged over a slightly 
smaller longitude range before unfolding. This may be an indication that the 
negative emissivities result in part from local fluctuations or sources, since 
the high energy data, being more confined to the plane (Mayer-Hasselwander et, 
al. 1982), is relatively less sensitive to local effects. 

IV. COSMIC RAY DISTRIBUTION 

The results presented in the previous section give distributions of the 
total observed gamma-ray emissivity in the Galaxy. This emission has two 
basic types of component: 1) diffuse emission from cosmic rays intei acting 
with gas or with a low energy radiation background, and 2) point source 
emission from pulsars, accreting compact objects, etc. The total point source 
contribution is uncertain, but estimates of the emission level from pulsars 
(Harding 1981, Harding and Stecker 1981, Salvati and Massaro 1982) are in the 
range 15% - 20% or less. We note that pulsars are the only true galactic 
gamma-ray point sources known, even though there exist 30 or so "point 
sources" with angular sizes of 1° - 2" in the COS-B catalog (Swanenburg et. 


al. 1981), many of which could be diffuse. At least 80% - 90% of the 
emissivlty at low latitudes, where the Compton component Is low (Kniffen and 
Fichtel 1981), can then be presumed to result from Interactions between cosmic 
rays and gas In the form of molecular and atomic hydrogen. In addition, the 
galactic point source emission may be distributed like the diffuse emission 
(Harding and Stecker 1981) so that Its effect does not distort the overall 
flux distribution. Therefore, Information on the distribution of gas In the 
Galaxy can be used In conjunction with the observed gamma-ray emissivlty to yield 
Information on the galactic cosmic ray distribution (e.g, Stecker and Jones 
1977). 

The quantity, the gamma-ray emissivlty per H-atom, Is derived from 

the observed gamma-ray volume emissivlty, total gas density, njoT* 9^® 

scale height, hg, by the following: 


QJr)h 

(r) = — fp 

"tot "g 


(5) 


The total gas density Is the sum of molecular, n^^, and atomic, n^, 
densities, n^-Qj = 2 nj^^ ■*” nm* Molecular hydrogen densities have been derived 
from galactic CO surveys which have recently been extended to Include southern 
galactic longitudes. Longitude-velocity data from these surveys can be 
unfolded using a galactic rotation curve to give CO radial emissivlty 
distributions, which can then be converted to molecular hydrogen densities. 

We use here the recent CO survey results of Robinson et. al. (1984), who give 
radial distributions of J (K km s"^ kpc"M In the ^^0 (1-0) transition for 
both the North and South, and scale heights as a function of galactic 
radius. The radial CO distribution shows the bulk of emission concentrated In 
a ring of radius '' 5 kpc, with the emission In the South distributed In a 
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broad peak between 2.5 and 8,5 kpc, In contrast to ttie sharp peak in the North 
at 5 kpc. Using their conversion factor, n^^^ « 0.12 J (CO^^) cm"^, we 
computed the molecular gas densities in 1 kpc radial bins, and then added a 
constant atomic hydrogen density of 0.35 cm“® for R > 4 kpc. A combined scale 
height was obtained for each radial bin as a weighted average of the molecular 
scale heights and a constant HI scale height of 250 pc. The values for Hq 
calculated in this way increase by at least a factor of two between R = 2 kpc 
and R = 10 kpc due to both an increasing CO scale height and an increasing 
fraction of atomic gas, which has a much larger scale height than the CO. 

Figure 5 shows the radial distribution of > 100 MeV ) derived using 
the SAS-2 and COS-B emissivities plotted in Fig. 1. If all the gamma-ray 
emissivity were from diffuse processes, then qY would be proportional to the 
density of cosmic rays. The emissivity per H-atom derived from both the SAS-2 
and COS-B -how evidence for an increase in the inner galaxy relative to 
local valits in both the North and the South. The COS-B gradient, however, is 
significantly smaller than the SAS-2 gradient due to smaller COS-B emissivities 
in the inner Galaxy. In fact the value of qy is less than the local value 
over a number of bins, which would require a decrease in cosmic ray density to 
compensate for a gradient in gas density steeper than the gamma-ray emissivity 
gradient. Figure 6 shows the distribution of qy using the COS-B emissivity 
for the case of a higher background subtraction (cf. Fig. 2), which increases 
the gamma-ray volume emissivity gradient and thus the implied cosmic ray 
gradient, bringing it into agreement (within the errors) with the SAS-2 result. 

We can also compare the gamma-ray emissivity per H-atom for the COS-B 
data in the two energy ranges E > 100 MeV and 0.3 < E < 5 GeV. Since diffuse 
processes involving cosmic-ray protons produce different gamma-ray spectre 
than those processes involving electrons, variations in the cosmic-ray 
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gradient as a function of energy would reveal differences In the distribution 
of cosmic-ray protom; and electrons as suggested by Bloemen, et al . (1984) for 
the outer Gala>y. In the 0.3 - 3 GeV energy range the distribution of qy 
plotted in Fig. 7 shows no statistically significant variation from the lower 
energy distribution plotted In Fig. 5, Indicating that, at least in the Inner 
Gala;y^, cosmic-ray protons and electrons probably have the same radial 
distribution. 


V. DISCl'SSION OF RESULTS 

Our results generally confirm the phonomenology of the 5 kpc “Great 
Galactic Ring" and cosmic ray gradient derived earlier using the northern 
hemisphere data (Stacker, 1975, Stacker 1977) and extend the ring 

phenomenology to the southern hemisphere region of the Milky Way. In 
addition, the southern galactic region gives evidence of more structure than 
the northern region does. This Is true for the CO data, the Y-ray data, and 
the radio data on the distribution of HI I regions and pulsars In the Inner 
galaxy* In particular, the southern region gives evidence of enhanced Y-ray 
emission at a radial distance of 7 to 8 kpc from the galactic center (noted 
previously by Stacker 1977 and Carav'ao and Paul 1979). Such emission Is more 
evident In the COS B data than the SAS-2 data. This emission region 
represents only a small fraction of the emission coming from the 5-kpc ring. 

An interpolation of the northern and southern emlssivlty data as given in Fig. 
3 leads most naturally to the Interpretation that the emission Is associated 
with a spiral arm feature of the Galaxy, the Carina arm. The overall 
emlssivlty pattern shown In Fig. 3 Is consistent with the spiral pattern 
suggested by Georgelln and Georgelln (1976) based on the distribution of HI I 
regions. The main Y-ray emission appears to He In a region Inside the spiral 
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pattern of the largest giant molecular clouds derived by Thaddeus and co- 
wo'kers (to be published). All of these Indications support the thesis that 
the galactic Y-ray emission Is associated with the youngest regions of star 
formation In the Gala>o' (Stecker 1976). 

The existence of the 5 kpc ring of molecular clouds and associated 
stellar activity now appears to be beyond question. Molecular rings have also 
been seen In CO observations of a number of spiral galaxies (Young and 
Scovllle 1982 and references therein.) Young and Scovllle (1982) attribute 
their origin to exhaustion of gas by star formation In the nuclear bulge of 
spiral galaxies. It Is Interesting to note that prominent ring structures on 
a galactic scale are also produced In numerical simulations of the evolution 
of barred spiral galaxies. Such rings form at the ends of the bar and may be 
theoretically associated with the corotation radius (Schwarz 1984). 

It Is thus apparent that unfoldings of the y-ray data have provided us 
with Important clues as to the nature and structure of the Gala>^y, Independent 
of models based on data In other wavelength ranges. It Is equally true, 
however, that data obtained from radio and Infrared surveys which reveal the 
overall distribution of gas and dust In the Galajo' can be used In conjunction 
with the Y-ray data to derive a more complete synoptic picture (Stecker 
1981). In particular, we find here that the > 100 MeV data reveal evidence 
of a cosm1c-ray gradient in the Inner Galaxy on both the northern and southern 
hemisphere sides of the galactic disk. We find some evidence of this also in 
the highest energy COS-8 data, which has the best angular resolution, and 
which also may be the only region of the spectrum clearly dominated by pi on 
decay Y-rays (Stecker 1977). These highest energy data have taken on 
additional Importance with the claim of Bloemen, et al. (1984) that there may 
be no gradient 1n the nucleonic component In the outer Galaxy (a claim which, 
If true, would pose serious theoretical problems). 


The cosmic-ray gradient derived for both the northern and southern 
galactic hemispheres also appears to peak in the region of the Great Galactic 
Ring where the supernova and pulsar distributions also have a maximum. This 
may most naturally argue for a galactic origin of the bulk of the nucleonic 
cosmic Y-rays, which have energies in the 1-10 GeV range. It has, of course, 
been realized for a long time that the cosmic ray electron component must be 
of galactic origin (Fazio, Stecker and Wright 1966). 

We therefore feel that the existence of a radial gradient in the galactic 
cosmic ray distribution is now on firm ground. In fact, there is evidence 
that (owing to increased stellar nucleosynthesis) the ratio of CO to 
abundance may increase toward the galactic center, perhaps by up to a factor 
of 2 in the 5-6 kpc ring (Shaver, et al. 1983). This would reduce estimates 
of the average volume density in the inner Galaxy and would actually 
increase the derived gradient in the cosmic-ray flux by a corresponding 
factor. Such as increase would not qualitatively change our conclusions and 
would strengthen the argument for a general galactic cosmic-ray gradient. 
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FIGURE CAPTIONS 


Fig. 1. Gamma-ray surface emissivity (volume emissivlty times scale height) 
as a function of galactocentric radius In units of Rg = 10 kpc, derived from 
the SAS-2 and COS-B data at energies greater than 100 MeV. The right hand 
scale shows the emissivltles relative to the local value. 

Fig. 2. Gamma-ray surface emissivlty distribution derived from the COS-B data 
at energies greater than 100 MeV with an additional background subtraction 
(see text) . 

Fig. 3. Schematic map of the Y-ray emissivlty In the Inner Galaxy looking down 
on the galactic plane, constructed from the radial distributions shown In Fig. 
1, omitting the galactic center region. The positions of the Sun and the 
galactic center are inllcated by the appropriate symbols. 

Fig. A. Gamma-ray surface emissivlty distribution derived from the COS-B data 
In the high energy (0.3 - 5 GeV) channel. 

F1g. 5. Radial distribution of Y-ray emissivlty per H atom above 100 MeV in 
the Northern and Southern halves of the Galaxy derived from the radial surface 
emissivlty distributions shown in Fig. 1 using the CO data of Robinson, et al 
(1984). 

Fig. 6. Emissivlty per H atom as In Fig. 5, but using an extra background 
subtraction to the COS-B data (see text). The results are compared with those 
obtained from the SAS-2 data. 
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Ffg, 7. Gamma-ray emisstvity per H atom derived from the COS-B radial surface 
emissivity distribution in the high energy channel (0.3 - 5 GeV) shown in Fig. 4. 
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